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Abstract 

We study the stability of the color-flavor locked (CFL) phase of dense quark matter with regard 
to the formation of a non-zero Goldstone boson current. We show that an instability appears in 
the vicinity of the point /i^ = A which marks the appearance of gapless fermion modes in the CFL 
phase. Here, /x^ = m^/(2/u) is the shift in chemical potential due to the strange quark mass and A 
is the gap in the chiral limit. We show that in the Goldstone boson current phase all components of 
the magnetic screening mass are real. In this work we do not take into account homogeneous kaon 
condensation. We study the effects of an instanton induced interaction of the magnitude required 
to suppress kaon condensation. 
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I. INTRODUCTION 



Exploring the phase diagram of dense baryonic matter is an area of significant theoretical 
and experimental activity. Using methods based on weak-coupling QCD it was shown that 
the groundstate of three flavor quark matter at very high baryon density is the color-flavor- 
locked (CFL) phase 

BBS. 

Finding the second, third, etc. most dense phase turns out 
be a much harder problem. If the quarks are degenerate and flavor symmetry is exact then 
the symmetries of the CFL state match the symmetries of (hyper) nuclear matter. This 
means that high density quark matter and low density nuclear matter might be continuously 
connected, without an intervening phase transition |^. 

In the real world the strange quark mass is not equal to the masses of the up and down 
quark and flavor symmetry is broken. At high baryon density the effect of the quark masses 
is governed by the shift Hq = m^/(2/i) of the Fermi energy due to the quark mass. This 
implies that flavor symmetry breaking becomes more important as the density is lowered. 
Clearly, the most important effect is due to the strange quark mass. There are two scales 
that are important in relation to fis- The flrst is the mass of the lightest strange Goldstone 
boson, the kaon, and the second is the gap A for fermionic excitations. When /i^ becomes 
equal to tuk the CFL phase undergoes a transition to a phase with kaon condensation P|. 
In the present work we wish to study possible phase transitions in the vicinity of the point 
fis = A. 

Alford, Kouvaris and Rajagopal observed that gapless fermions appear in the spectrum 
if /is > A and that the CFL state is energetically stable compared to other homogeneous 
phases as long as fis < 2 A. They termed the phase in the regime A < fig < 2 A the gapless 
CFL (gCFL) phase 0, Q| • The problem is that gapless fermion modes in weakly coupled 
superfluids tend to cause instabilities in current- current correlation functions 0, 0] and these 



instabilities have indeed been found in the gCFL phase 



Id, 
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We have suggested that 

the instability is resolved by the formation of a non-zero Goldstone boson current jl^l . Q| . 
We will refer to this phase as the current-CFL (curCFL) phase. 

In the present work we wish to show that the spacelike screening (Meissner) masses in the 
curCFL phase are real, and that the instability is indeed resolved. We also improve on the 

ways. We allow the CFL gaps Ai^2,3 to be different and we 
properly implement electric charge neutrality. On the other hand, we make one important 
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simplification as compared to our earlier work and ignore homogeneous Goldstone boson 
condensates. This assumption simplifies the calculation of the dispersion relations and allows 
a more detailed comparison with previous work on the gCFL phase. We will comment on 



the effect of kaon condensation in the cone 
the LOFF phase in three flavor QCD 



usions. We a 
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currents in two flavor QCD can be found in 
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2G 



21 



so compare our results to studies of 
. Recent work on Goldstone boson 



Jl9LR( 
22, 2^. 



II. EFFECTIVE LAGRANGIAN AND DISPERSION LAWS 



We consider an effective lagrangian t 



with background gauge flelds 
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m 



lat describes the interaction of gapped fermions 



C = Tt [xii^v -d-fi"^- A,Q)xl) + IV (xi{^v ■ 9 - /i« - A,Q)xr 



(1) 



a,b,i,j,k 



Here, xT r. left /right handed fermions with color index a and flavor index i, are 
SU{'i)c color gauge flelds, and fi^ = MMV(2/i), /t^ = M'^M/{2jd) are effective chemical 
potentials induced by the quark mass matrix M. The matrix Q = diag(|, — |, — |) is the 
quark charge matrix and A^ is an electro-static potential. The flelds X, Y determine the 
flavor orientation of the left and right handed gap terms and transform as X — LXC^, 
Y RYC^ under (L, R) G 5f/(3)L x SU{'i)R and C G 5f/(3)c, and [k = 1, 2, 3) are 
the CFL gap parameters. The gauge fleld A^ appears transposed because of the transforma- 
tion properties of X, Y. From the lagrangian equ. ([Q) we can read off the Nambu-Gor'kov 
propagator (left handed part) 




(2) 



where p = v-p—fi, with the Fermi velocity v. The components of the propagator are matrices 
in color-flavor space. We use a basis spanned by the Gell-Mann matrices {A = 1, ... ,8) 
and Ao = \/fl- In this basis 



A 



AB 



lAabeijaersb>^{l.>^fs, ^ab = diag(Ai, A2, A; 



(3) 
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We will assume that X = Y = 1 which excludes the possibility of kaon condensation. This 
assumption significantly simplifies the calculation of the fermion dispersion relations and 
the current correlation functions. It is possible to suppress kaon condensation by including 



a large instanton induced interaction 



26| . We will discuss instanton effects in Sect. IVll The 



gauge field vertex in the phase with X = y = 1 is 



'^v,AB = |Tr IXAifi^ + AeQ)XB + XaXbv ■ 



(4) 



We define 



A 



or 



-/i^ - A,Q + A3A3 + AgA, 



8; 



(5) 



where /I3, A^ are the diagonal components of the color gauge field relative to the leading 
order solution of the color neutrality condition. We wish to study the possibility of forming 
a Goldstone boson current. Gauge invariance implies that the free energy only depends on 
the combinations Jl = X'^{S/ — iA^)X and Jr = ^"'^(V — iA'^)Y. We will restrict ourselves 
to diagonal currents Jl.r- Within our approximations the free energies of the vector and 
axial-vector currents Jl = ^Jr are degenerate. We will consider the pure vector current 
Jl = Jr = A^ with 

^^^ = i/(^3 + ;^As-^Ao). (6) 

This ansatz has the feature that it does not shift the energy of electrically charged fermion 
modes. We will see that an ansatz of this type is favored if electric neutrality is enforced. 
The electric charge due to a gapless charged fermion is 6Q ~ /i^j. The density of electrons is 
Ue ~ /ig and electric neutrality gives /ig ~ (i^'^jY^^ ^ J- This means that any solution with 
fif, J A has to have a very small density of gapless charged fermions. We will study 
a more general ansatz for the current in Sec. |3 and show that the energetically preferred 
solution is indeed close to the ansatz equ. (jH))- From the propagator (0) we obtain the 
following dispersion laws, 



€2 



A3 

2 



VSA^ 



V ] 



+ 



A2 



h H h 

2 2 2 



' _^Aa Aa^ 



A2 



^3 



+ A^- fXs + 
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^8 V-f 
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-A3 + + 



As V-j\ 



p2 + A? 



p2 + Af + . . . , 



p2 + 4Af + . . . . 



(7) 



The dots in ey^g.g denote terms of the order a" {n > 1) with a G {^3, Ag, t?- J, A2 — Ai, A3 — 
Ai}. We observe that without a current ei and €3 become gapless at fig = A. We also note 
that the current shifts the energy of the electrically neutral modes ei 2 and the momenta of 
the charged modes e3_6. 

III. FREE ENERGY 



In order to determine the groundstate we need to supplement the effective lagrangian 
with a potential for the order parameter Aj^.. In QCD this potential is determined by gauge 
field dynamics. For the present purpose we are not interested in the details of this interaction 
and we will model the effective potential as a quadratic function of Afc. The free energy is 
given by 



G 



(A2 + A^ + A 



1 



2„2 



87r2 



47r2 
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ldp[dtY^{\e,\-\p\), (8) 

i=l 



with t = cos 6. Here, G fixes the magnitude of the gap in the chiral limit and Gi takes into 
account instanton effects. The structure of the instanton term is determined by vacuum ex- 
pectation value of the 't Hooft interaction. The value of Gi can be computed in perturbative 
QCD Q but in the present work we will treat it as a free parameter. We first study the case 
that there are no instanton effects, = 0. The integral in equ. (jHl) is quite complicated. 
We assume that the external fields and currents are small and expand in a set of parameters 
b e {is/Ai, ig/Ai, j/Ai, (A2 - Ai)/Ai, (A3 - Ai)/Ai, (/i, - Ai)/Ai}. We will study the 
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convergence properties of this expansion below. To second order in b we obtain 

n c i(A? + A^ + Al) + ^[-54(A? + A^ + A^)log(|-) 

+ A^(15 -341og2) + 2(A^ + A^)(21 - 171og2) 

- I6A1 ( A2 + A3) (6 - log 2) + 4A2 A3 (3 + 4 log 2) 

2/^2 



+ ^ (219 - 32 log 2) + ^ {Al + Al){-3 + 2 log 2) 



6487r2 

I ' 2 2^3' 



As A3 ^ VSAs 



A, 



2 2^3' 2 

A, 



Ae + IJ,s,J,0 



+ HlA3,-^,0,j,0 



where the function H is given by 



H{A,ai,a2,ji,j2) = 



127r2 



6 AM 1 + 2 log 



/2^ 
V A 



+ 12a2+j2 



+ 



2^A/A 
157r2_72 



(2a2 - 2A - j2)'/'e(2a2 - 2A - J2) 
(2a2 - 2A + j2)'/'e(2a2 - 2A + J2) 



The gap equations and neutrality equations can be summarized as E'j = 0, where 

dfl , 



E4 



a A/ 
dfl 



OA,' 



z = 1,2,3) 
E5 — 7n~5 Eq 



OAs' dAs 

Prom the hnear combination E4 + + we obtain the relation 

is = -VSAs 

Including the {—A^) term in fi, we find from E4 (for n ^ A^) 

Ae = max(0, - A2 - A3). 



(9) 



(10) 

(11) 
(12) 

(13) 
(14) 



We assume — A2 — A3 < (and therefore Ae — 0), which will turn out to be consistent 
for the values of fig under consideration. Then we find A3 = A2. Furthermore, E2 is then a 
linear equation in A2, with the solution 

54y \ 



A2 = Al 1 



51 + 54y-41og2y ' 



(15) 
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FIG. 1: Left panel: O/(/x2A(0)2) as a function of j/A(0) for Ai = A2 = A3 = A(0), Ae = 
ig = ig = at = 0.9683A(0) {0{b'^) approximation). Right panel: function 
of j/A(0) with gap equations and neutrality conditions taken into account. The different curves 
correspond to different orders in the expansion in b: second order (black), third order (green), 
fourth order (blue), fifth order (red). All curves are shown for the respective values of i^s^crit- 

where y = log(Ai/A(0)) and A(0) is the value of the gap in the chiral limit. The linear 
combination Ei — E4 — 2E^ gives a linear equation for ^3, with the solution 



where s± = 2^s — 2Ai — AA^ ± j. We solve equ. (fT7|) numerically for Ai for given values of 
fig and J. 

We observe that Q/{fi'^A{Of) is a function of 6/A(0) and Ai/A(0). This means that 
in order to study the effective potential as a function of dimensionless variables we do not 
have to specify the value of fi and G. We first study the dependence of the thermodynamic 
potential Q on the current j for different fis- We begin with the simplified case Ai = 
A2 = A3 = A(0), Ae = A3 = As = 0. We find that f2(j) develops a nontrivial minimum 
for /is > 0.9683A(0), as shown in the left panel of Fig. [T] If we take into account the gap 
equations and neutrality conditions, a nontrivial minimum appears for a slightly larger value 
of fis, namely /i^ > fis,crit = 0.9916A(0), as shown in the right panel of Fig. [T] (black curve). 
We note that the expectation value of j is considerably smaller than in the simplified case. 

We have also studied the importance of higher order terms in b in the free energy. These 
terms are not included in equ. ()9I1U|) but they can easily be generated using a symbolic 




(16) 



There remains one equation (proportional to E4 + 2E5), 




(17) 
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FIG. 2: Optimum current j/A(0) as a function of ;Us/A(0). We compare results at different orders 
in the small parameters b (curves labeled as in Fig. 

algebra package such as Mathematica or Maple. In Fig. ^we compare the 0(6^) result to 
the 0(6") approximation with n = 3,4, 5. At the minimum of fl we find A3 = A2 and Ae = 
as before. We also see that while the second order approximation gives results that are only 
qualitatively correct, the differences between the results of the higher order approximations 
are very small. 

We may also determine the minimum of the free energy for fig greater than the critical 
value. Fig.Elshows the expectation value of j as a function of fig. Again, we compare results 
at different order in the small parameter b. We can see that the differences between the 
results of the higher order approximations are very small even for (somewhat) higher values 
of fis- The highest approximation we discuss is the fifth order approximation. In that case 
we find that the critical values of fig is given by fis,crit = 0.9919A(0). The expectation value 
of the J is equal to 0.1541A(0) at this point. 

In Fig. iniwe show the gauge potential A3 and the flavor asymmetries in the gaps, Ai — 
A(0), A2 — A(0), as a function of fig calculated to 0(6^). We remark that ei becomes a 
gapless mode for /x^ > /^s,crit, while 63 (the almost quadratic mode of Ref. Q) has at least a 
very small gap, greater than 0.04A(0), for all values of fig. 

IV. MEISSNER MASSES 

In this section we study the stability of the current-current correlation functions. The 
screening (Meissner) masses are given by 
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FIG. 3: Gauge field A3/A(0) (dotted) and flavor symmetry breaking in the gap parameters (Ai — 
A(0))/A(0) (continuous) and (A2 — A(0))/A(0) (dashed) as functions of ^s/A(0), computed to 
0{b^). 



X Tr 



g+{p)v:g+{p + k)vi + g-{p)v:g-{p + k)Vi 



+ s+(p)k:s-(p + k)vi + h-(p)v:h+(p + kwi 



where we have defined the vertices 



with -y* = —v^. We begin by studying the screening masses for sniaU /i^ < /i^ 



,crit- 



;i8) 



(19) 



In this 



case the Meissner mass matrix is diagonal, and we have m 
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= m\.j^j, with 
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1 /a^ 
— - arctan — 



1 



a 



■arctanh 
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(20) 

(21) 



with a = /is/A(0), /? = V4 — and 7 = V9 — lOa^ + q;4_ We note that the Meissner 
masses of the modes i = 1,2,3,8 are independent of fis below fis,crit- These modes only 
couple to pairs of quasi-particles that are subject to the same shift in energy. In our ap- 
proximation, in which we neglect the effect of fig on the density of states, this implies that 
the /is dependence can be eliminated by a simple shift in the loop energy. This is not the 
case for the off-diagonal gluon modes i = 4,5, 6, 7. We find that the screening masses rriMn 
{i = 4, 5, 6, 7) drop with fig, but they remain real and non-zero for /i^ < ^s,crit jl^ - 

In the presence of a finite current j we may decompose the Meissner masses into a 
longitudinal and a transverse component. 
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FIG. 4: Meissner masses 



squared m^^^^^^J^g ^ ), mj^,^(2)||/(5 ^ )' "^Mii±/(5 M ) 



and 



m 



M\\\\ 



\\\l{.9^^) as functions of /is/A(0), computed to 0(6^). 



First let us evaluate the Meissner masses through order 6°, which corresponds to the 0(6^)- 
approximation in the free energy. We will see that the interesting terms are those of order 

A chromomagnetic instability could occur for the Meissner masses with color indices 
1,2,3 or 8 Q]. In the 3-8-sector there are contributions where both quark propagators in 
the gluon self energy diagram are gapless. We find that the mass matrix in the 3-8-sector is 
not diagonal, the mixing angle being equal to |. We obtain the following eigenvalues of the 
3-8-mass matrix within our approximations, 



M(l)±,|| 



2 2 21 - 81og2 
108.2 ' 

2 2f 21-81og2 
1087r2 
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^1 



457r2j3 



0(s+)sf (2s+-5j) 
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M{2)|1 



2 2f 21-81og2 
108vr2 



907r2/ 

e(s_)^(15/ + 20js„ + 8s^)]|. 



[e(s+)y57(15/-20js+ + 8s^ 



(23) 



(24) 



(25) 



These quantities are positive for /i^ = fis,crit when j is chosen such that Q is minimal. We 
remark that the longitudinal Meissner mass squared fn\,j^2)\\ is very similar to the second 
(partial) derivative of Q with respect to j (only the G-independent terms are different because 
the current contains a component proportional to Aq.) 

For m\j^^{= m\j22) there are contributions where one quark propagator in the gluon self 
energy diagram is gapless and the other one contains the the would-be quadratic mode. We 
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find 
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e(s+) v^(8s; + 15r^ + 25rs, 



(j ^ -J, s+ ^ s- 



(27) 



where r = 2(A2 — /is + A3) > 0. Again we find that "ii/nx and positive. 

The squared Meissner masses m|^44 to m%jjY equal at 0(6°). At the above value of 
/is we find 

2 9-y37r-31og2 
mj\^44 = ^^^^^ , (28) 

which is again positive. 



As above we may include higher orders in b in the calculation of the Meissner masses. 
Fig. m shows the "dangerous" components of the Meissner masses as functions of fis at 
0{b^). We observe that all masses are real. As above the differences between the results of 
the higher order approximations are very small. The numerical values of the "dangerous" 
components of the Meissner mass just above the phase transition are shown in the first line 
of Table m 



V. AN IMPROVED ANSATZ FOR THE CURRENT 

We have studied the effect of using a more general ansatz for the current. Consider 

= (ciAs + cs^Ag - ca^Ao) (29) 

The ansatz given in equ. © corresponds to ci = C2 = C3 = 1. We now minimize fl with 
respect to j and the coefficients Cj. We can fix ci such that the explicit term in equ. (jS)) 
remains unchanged as C2 and C3 are varied, which amounts to setting ci = y^(9 — 2c2 — c^/G. 
At 0{b^) we find Hs,crit = 0.9918A(0), and 

ci = 0.980, C2 = 1.000, cs = 1.114. (30) 
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"^mii||Mm(0) 


Ci = l 


0.356 


0.056 


0.707 


0.174 


Cl,3 / 1 


0.358 


0.061 


0.696 


0.142 



TABLE I: Dangerous components of the Meissner masses at ^u^ = fJ,s,crit in units of the screening 
mass at = (third order approximation). First hne: ci^2,3 = 1 in equ. second hne: Cj 

according to equ. (|5n|) . 

We observe that the exact minimum is indeed quite close to our simple ansatz. The "dan- 
gerous" components of the Meissner mass are not significantly affected in the vicinity of the 
phase transition. Their values just above the phase transition are shown in the second line 
of Table m 



VI. INSTANTON CONTRIBUTION 



In this section we consider the role of an instanton induced interaction. The value of 
Gi can be computed in perturbative QCD 2^ but the result is very sensitive to the QCD 
scale parameter and to higher order perturbative corrections. In the present work we have 
neglected kaon condensation and we will therefore explore the effect of an instanton term 

n 

of the magnitude required to suppress kaon condensation. We have previously argued |26|| 
that an instanton terms of this size is problematic from the point of view of instanton 
phenomenology at zero baryon density, but we will ignore this issue here. 
In order to suppress kaon condensation, has to satisfy 

1 f2u3/2 



where is the pion decay constant in the CFL phase. We have used_ 
together with the perturbative result for the pion decay constant 



/G, = 0.67/,^A(0) 



2 _ 21-81og(2) / ^2 



In this case the effective potential develops a non-trivial minimum for /ig > ^s,crit = 
0.964A(0) (using n to 0(6^) and d = 1). At this point we find A3 = 1.21A(0), such that the 
condition in equ. is just met. The magnitude of the current at the onset is j = 0.17A(0) 
and the shape of the effective potential is not strongly modified by instanton effects. We also 
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FIG. 5: Aj/A(0) as functions of /Us/A(0) with the instanton contribution taken into account 
{l/Gi = 0.67/2A(0)"1/2)^ computed to 0{h^). Below the phase transition Ai and A2 are identical. 

find that all Meissner masses are real. The main effect of the instanton term is a non-zero 
splitting between A3 and Ai 2 below fis,crit, as shown in Fig. El This effect leads to the 
small reduction in the critical value of the effective chemical potential, Hs,crit = 0.964A(0) 
compared to fis,crit = 0.992A(0) without the instanton term. We also note that because of 
the larger splitting between the gaps our expansion converges more slowly. 

VII. CONCLUSIONS 

We have shown that near the point at which gapless fermion modes appear in the spec- 
trum the CFL phase becomes unstable with respect to the formation of a Goldstone bo- 
son current. We have computed the magnetic screening masses in the Goldstone current 
phase and shown that all masses are real. Our calculation is based on an effective la- 
grangian of fermions coupled to a mean-field gap term and background gauge potentials. 
The CFL pairing energy is of order /i^A^. We have performed an expansion in powers of 
{{fig — A)/A, y4/A, j/A, 5A/A) and included all orders in this expansion until numerical 
convergence was obtained. We have neglected all contributions that are suppressed by extra 
factors of A/fi. 

In this paper we have not included the effect of a homogeneous kaon condensate. Kaon 
condensation is likely to play an important role since the kaon condensate modifies the 
spectrum of fermions in such a way that the lightest mode is charged. We have considered 
the possibility that the condensate is suppressed by a large instanton term. We find that 
the instanton term leads to a splitting between the different CFL gaps for fig < fis,crit but it 
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does not qualitatively change the effective potential for the Goldstone boson current. 

The phase studied in this paper is equivalent to a single plane wave LOFF state. The 
constant gauge field A can be removed by a gauge transformation which leads to the following 
gap matrix (for q = 1), 

A,, = diag(Aie-*^^-^ A^, Aa)^^. (33) 
This gap matrix is different from the ansatz Aab = diag(0, A2, A2)afee~*^'-'* considered by 



Casalbuoni et al. in 



13] . One important difference is that we consider the CFL phase near the 



onset of the Goldstone boson current instability. In that regime all fermions remain paired 
and the free energy can be computed as an expansion in (j/A). Casalbuoni et al. consider 
the CFL phase near the maximum value of fig at which pairing between different flavors is 
possible and expand the free energy in (A/j). Clearly, it is important to understand how to 
interpolate between these two extremes. Once a kaon condensate develops the current is no 
longer diagonal, and it is also not equivalent to a pure gauge field. 
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